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Abstract 

Confinement in QCD results from special properties of vacuum 
fluctuations of gluon fields. There are two numerically different scales, 
characterizing nonperturbative QCD vacuum dynamics: "small" one, 
corresponding to gluon condensate, critical temperature etc, which 
is about 0.1-0.3 GeV, and a "large" one, given by inverse confining 
string width, glueball and gluelump masses etc, which is about 1.5- 
2.5 GeV. We discuss this mismatch in a picture where confinement is 
ensured by quadratic colorelectric field correlators of the special type, 
denoted by D^{x). The latter together with four other quadratic 
correlators, , Df , and D^^ support most of perturbative and 
nonperturbative dynamics, while contribution of higher order terms 
is shown to be small. These correlators, on the other hand, can be 
calculated via gluelump Green's function, whose dynamics is defined 
by the correlators themselves. In this way one obtains a self-consistent 
scheme, where string tension can be expressed in terms of Aqcd- Thus 
the dynamics of vacuum gluon field fluctuations supports mean-field- 
type mechanism leading to confinement at temperatures T <Tc. 



1 Introduction 



Confinement of color is tlie most important property of quantum chromo- 
dynamics (QCD), ensuring stability of matter in the Universe. Attempts to 
understand physical mechanism of confinement are incessant since advent of 
constituent quark model and QCD, for reviews see e.g. [U El El H] . Among 
many suggestions one can distinguish three major approaches: 

• confinement is due to classical field lumps like instantons or dyons 

• confinement can be understood as a kind of Abelian-like phenomenon 
according to seminal suggestion by 't Hooft and Mandelstam [5] of 
dual Meissner scenario. This approach has gained much popularity in 
the lattice community and studies of various projected objects such as 
Abelian monopoles and center vortices are going on (see, e.g. [6]). 

• confinement results from properties of quantum stochastic ensemble of 
nonperturbative [np) fields filling QCD vacuum. The development of 
this approach in a systematic way started in 1987 [7] (see, e.g. [8] 
as example of earlier investigations concerning stochasticity of QCD 
vacuum). The nontrivial structure of np vacuum can be described by 
a set of nonlocal gauge-invariant field strength correlators (FC). We 
discuss this scenario in the present paper. 

QCD sum rules P, (TUj were suggested as an independent approach to 
npQCD dynamics, not addressing directly confinement mechanism. The key 
role is played by gluon condensate G2, which is defined as np average of the 
following type: 

In the sum rule framework it is a universal quantity characterizing QCD 
vacuum as it is, while it enters power expansion of current-current correlators 
with channel-dependent coefficients. This is the cornerstone of QCD sum 
rules ideology. 

It is worth stressing from the very beginning that G2 and other conden- 
sates are usually considered as finite physical quantities. Naively in pertur- 
bation theory one would get G2 ~ a~'^, where a is a space-time ultraviolet 
cutoff (e.g. lattice spacing). It is always assumed that this "hard" contribu- 
tion is somehow subtracted from ([T]) and the remaining finite quantity results 
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from "soft" np fields, in some analogy with Casimir effect where modifica- 
tion of the vacuum by boundaries of typical size L yields nonzero shift of 
energy- momentum tensor by the amount of order L~^. In QCD the role of 
L is played by dynamical scale ^qcd- 

The idea of np gluon condensate has proved to be very fruitful. How- 
ever the relation of G2 to confinement is rather tricky. Let us stress that 
since there can be no local gauge-invariant order parameter for confinement- 
deconfinement transition, G2 is not an order parameter and, in particular, 
neither perturbative nor np contributions to this quantity vanish in either 
phaseJll On the other hand one can show (see details in [3]) that the scale 
of deconfinement temperature is set by the condensate, or, more precisely, 
by its electric part: ~ Let us remind the original estimate for the 

condensate given by G2 = 0.012 GeV^, with large uncertainties. One 
clearly sees some tension between this "small" scale and a "large" mass scale 
in QCD given by, e.g. the mass of the lightest 0"*"+ glueball, which is about 
L5 GeV. 

The question about the origin of this mismatch turned out to be rather 
deep. It was suggested in [7], see also [HI |T2] that there is another impor- 
tant dimensionfull parameter characterizing np dynamics of vacuum fields: 
correlation length A (also denoted as Tg in some papers). The crucial feature 
of stochastic picture found in [13] is that the lowest, quadratic, nonlocal FC, 
{F^^{x)Fxa{y)) describes all np dynamics with very good accuracy. It was 
also shown that a simple exponential form of quadratic correlators found on 
the lattice pHj allows to calculate all properties of lowest mesons, glueballs, 
hybrids and baryons, including Regge trajectories, lepton widths etc, see re- 
views [21 [12] and references therein. Since the correlation length A entering 
these exponents, is small, potential relativistic quantum-mechanical picture 
is applicable and all QCD spectrum is defined mostly by string tension a 
(which is an integral characteristic of the nonlocal correlator, see below) and 
not by its exact profile. This is discussed in details in the next section. 

However to establish the confinement mechanism unambiguously, one 
should be able to calculate vacuum field distributions, i.e. field correlators, 
self-consistently. In the long run it means that one is to demonstrate that 
it is essential property of QCD vacuum fields ensemble to be characterized 
by correlators, which support confinement for temperatures T below some 
critical value Tc, and deconfinement at T > Tc. 

^Let us repeat that following ^ we define G2 in ([T]) as purely np object. 
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Attempts to achieve this goal have been undertaken in [TJ], however the 
resulting chain of equations is too complicated to use in practice. 

Another step in this direction is done in [16], where FC are calculated 
via gluelump Green's functions and self-consistency of this procedure was 
demonstrated for the first time. These results were further studied and con- 
firmed in [T7]. 

The main aim of this paper is to present this set of equations as a self- 
consistent mechanism of confinement and to clarify qualitative details of the 
FC - gluelump connection. In particular we demonstrate how the equivalence 
of colormagnetic and colorelectric FC for T = ensures the Gromes rela- 
tion [18J. We also show, that self-consistency condition for FC as gluelump 
Green's function allows to connect the mass scale to a^, and in this way to 
express Aqcd via string tension a. 

Lattice computations play important role as a source of independent 
knowledge about vacuum field distributions. Recently a consistency check of 
this picture was done on the lattice [19] by measurement of spin-dependent 
potentials. The resulting FC are calculated and compared with gluelump 
predictions in [20] demonstrating good agreement with analytic results, in 
particular small vacuum correlation length A ~ 0.1 Fm is shown to corre- 
spond to large gluelump mass Mq ~ 2 GeV. 

The paper is organized as follows. In the next section we remind the 
basic expressions for the Green's functions of qq and gg systems in terms of 
Wegner- Wilson loops and qualitative picture of FC dynamics is discussed. 
Section 3 is devoted to the expressions of spin-dependent potentials in terms 
of FC. We also shortly discuss the lattice computations of FC. In section 4 
FC are expressed in terms of gluelump Green's functions, while in section 
5 the self-consistency of resulting relations is studied. Our conclusions are 
presented in section 6. 

2 Wegner- Wilson loop, Field Correlators and 
Green's Functions 

Since our aim is to study confinement for quarks (both light and heavy) 
and also for massless gluons, we start with the most general Green's func- 
tions for these objects in a proper physical background using Fock-Feynman- 
Schwinger representation. The formalism is built in such a way, that gauge 
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invariance is manifest at all steps. A reader familiar with this technique can 
skip this section and go directly to the section 3. 

First, we consider the case of mesons made of quarks while gluon bound 
states are discussed below. For quarks, one forms initial and final state 
operators 

^ in, out 

where ip'^yip are quark operators, F is a product of Dirac matrices, i.e. 1, 
7/^,75, (7^75), and = Pexp{ig A^dz^) is parallel transporter 

(also known as Schwinger line or phase factor). The meson Green's function 
can be written (in quenched case) as 

Gqq{x,y\x\y') = {^^^t{x' ,y')^in{x,y))A,q = 

= (Ft$(x', y')Sqiy', y)F<|.(y, x)S,(x, x'))a (3) 
and the quark Green's function Sq is given in Euclidean space-time (see [21] 



and references therein): 

S,{x,y) = {m + D)-' = 

poo px 

= {m-D) ds{Dz)^yexp{-K)PAexp{ig A^dz^)P„{x,y, s), (4) 

Jo Jy 

where K is kinetic energy term, 

m is the pole mass of quark and quark trajectories z^{t) with the end points 
X and y are being integrated over in {Dz)xy 

The factor Pa{x, y, s) in (jl]) is generated by the quark spin (color-magnetic 
moment) and is equal to 



P„{x,y,s) = Ppexp 



9 / o-^^Ff,i,{z{T))dT 
. Jo 



(6) 



where a^,, = ^(7m7i'-7!/7m)' Pf{Pa) in ([6]), (Hj) are, respectively, ordering 
operators of matrices F^i,{A^) along the path z^ij). 
With ([3]), (HD one easily gets: 

Gqq{x,y]A) = 
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poo /"OO 

= ds ds\Dz),y{Dz%ye-''~'''Ti{T{m-b)W„{x,y)T\m! -b')), 
Jo Jo 

(7) 

where Tr means trace operation both in Dirac and color indices, while 

W^{x,y) = PAexp (ig A^dzAp^{x,y,s)Pl{x,y,s'). (8) 

V Jcix,y) / 

In ([8]) the closed contour C{x, y) is formed by the trajectories of quark z^{t) 
and antiquark z'^{t'), and the ordering Pa and Pp in P^, P^ is universal, i.e., 
^a{x,y) is the Wegner- Wilson loop (W-loop) for spinor particle. 

The factors (m — D) and {m' — D') in ([7]) need a special treatment when 
being averaged over fields. As shown in Appendix 1 of [22], one can use 
simple replacement, 

A - ^ f dZa\ 

m - D ^ m - ip, = 2 J " 

The representation Eq. ([7]) is exact in the limit Nc ^ oo, when internal 
quark loops can be neglected and is a functional of gluonic fields A^, F^^, 
which contains both perturbative and np contributions, not specified at this 
level. 

The next step is averaging over gluon fields, which yields the physical 
Green's function Gqq-. 

Ggq{x,y) = {Ggg{x,y-A))A. (10) 

The averaging is done with the usual Euclidean weight exp(— Action//?,), con- 
taining all necessary gauge-fixing and ghost terms. 

To proceed it is convenient to use nonabelian Stokes theorem (see [2] for 
references and discussion) for the first factor on the r.h.s. in ([8]) (correspond- 
ing to the W-loop for scalar particle) and to rewrite it as an integral over the 
surface 5" spanned by the contour C = C{x,y): 



{W{C)) = (^TrPexp \^tg J A^{z)dz^j j = (^TrPexp ^ig J F{u)ds{u 

oo „ oo 

= TrV. exp Y^itg}^ / {{F{u,) ..F{ur,)))ds,..ds,, = exp ^^"H*^] 

n=2 n=2 



5 



Here F{ui)dsi = ^{xQ,Ui)F^^{ui)t'^^{ui, Xojds ^y{ui) and Wj, Xq are the points 
on the surface S bound by the contour C = dS. The double brackets ((...)) 
stay for irreducible correlators proportional to the unit matrix in the colour 
space (and therefore only spacial ordering Vx enters ffTTj) ). Since ffTTl) is 
gauge-invariant, one can make use of generalized contour gauge |23], which 
is defined by the condition $(xo,Ui) = 1. Notice that throughout the paper 
we normalize trace over color indices as Trl = 1 in any given representation. 

Since (11 II) is an identity, the r.h.s. does not depend on the choice of the 
surface, which is integrated over in ds^y{u). On the other hand it is clear 
that each irreducible n-point correlator {{F...F)) integrated over S (these are 
the functions A*^")[S'] in (ITT]) ) depends, in general case, on the choice one has 
made for S. Therefore it is natural to ask the following question: is there 
any hierarchy of A'-"''[S'] as functions of n for a given surface S? 

To get the physical idea behind this question, let us take the limit of small 
contour C. In this case one has for the np part of the W-loop in fundamental 
representation: 

-log(iy(C))„p~G2-5'2 (12) 

where 5* is the minimal area inside the loop C. The np short-distance dy- 
namics is governed by the vacuum gluon condensate G2, see ([T]). Higher 
condensates and other possible vacuum averages of local operators, in line 
with the Wilson expansion, have been introduced and phenomenologically 
estimated in [10], for a review see p^ . 

Suppose now that the size of W-loop is not small, i.e. it is larger that 
some typical dynamical scale A = (9(1 GeV) to be specified below. Let us 
now ask the following question: if one still wants to expand the loop formally 
in terms of local condensates, what would the structure of such expansion 
be? It is easy to see from (fTT!) that there are two subseries in this expansion. 
The first one is just an expansion in n, i.e. in the number of field strength 
operators. It is just 

{W{C)) = 1 + A(2)[5] + A(3)[5] + ^(A(2)[5])2 + ... (13) 

The second series is an expansion of each A^"-* [S] in terms of local conden- 
sates, i.e. expansion in powers of derivatives. Taking for simplicity the lowest 
n = 2 term, this expansion reads 

A(2)[5]= / dsx [ dSy{{F{x)F{y))) = 
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= [ ds, [ dsy[{{Fix)Fix))) + ix-y),{{F{x)D,Fix))) + 

+ (x - y),{x - yU{F{x)D,D^F{x)))/2 + ...] (14) 

It is worth mentioning (but usually skipped in QCD sum rules analysis) 
that the latter series contain ambiguity related to the choice of contours in 
parallel transporters. The expression f[T^ is written for the simplest choice 
of straight contour connecting the points x and y (and not x and Xq or y 
and xq). This is the choice we adopt in what follows. In general case each 
condensate in f|T^ is multiplied by some contour-dependent function of x, y 
and xq. It is legitimate to ask whether this approximation, i.e. replacement 
of 

(Tr$(a;o, x)F^^{x)^x, a;o)$(xo, y)Fp^{y)<!>{y, xq)) 

by 

(TVF^,(x)$(x,i/)F,.(l/)$(l/,x)) 

affects our final results. Since the only difference between the above two 
expressions is profiles of the transporters, this is the question about contour 
dependence of FCs. We will not discuss this question in details in the present 
paper and refer an interested reader to [7j. Here we only mention that the 
result (weak dependence on transporters' choice) is very natural in Abelian 
dominance picture since for Abelian fields the transporters exactly cancel. 

Let us stress that both expansions (fT3|l and (fT4l) are formal since the loop 
is assumed to be large. Moreover, from dimensional point of view, the terms 
of these two expansions mix with each other. For example, one has nonzero 
v.e.v. of two operators of dimension eight: 

A^''>[S]^{F,,D'F,,) ; A^''>[S]^{F,^F,,D'F,,) (15) 

with no a priori reasoiJl to drop any of them. Moreover, even if one assumes 
that all these condensates of high orders can be self-consistently defined anal- 
ogously to G2, it remains to be seen whether the whole series converges or 
not. 

Here the physical picture behind PC method comes into play. It is as- 
sumed (and confirmed a posteriori in several independent ways) that if the 
surface S is chosen as the minimal one, the dominant contribution to (ITT!) 

^One could have such reasons if there is some special kinematics in the problem, but 
this is not the case. 
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results from v.e.v.s of the operators FD^F and this subseries can be summed 
up as (fT^ . In this language nonlocal two-point FC can be understood as a 
representation for the sum of infinite sequence of local terms f lT^ . 

Physically the dominance of (FD'^F) terms (and hence of two-point FC) 
corresponds to the fact that the correlation length is small yielding a small ex- 
pansion parameter ^ = FX^ ^ 1 (F is average modulus of np vacuum fields) 
or, in other words, typical inverse correlation length characterizing en- 
semble of QCD vacuum fields is parametrically larger than fields themselves. 
The dimensionless parameter ^ in terms of condensate is given by (6*2^^)^''^, 
and it is indeed small: of the order of 0.05 according to lattice estimates. 
One can say that often repeated statement " QCD vacuum is filled by strong 
and strongly fluctuating chromoelectric and chromomagnetic fields" is only 
partly correct - namely, in reality the fields are not so strong as compared to 
A~^. Technically one can say that we sum up the leading subseries in ffTTl) 
in the same spirit as one does in covariant perturbation theory [25] for weak 
but strongly varying potentials. 

Thus, for understanding of confining properties of QCD vacuum dynamics 
not only scale of np fields given by G2 is important but also another quantity 
- the vacuum correlation length A, which defines the nonlocality of gluonic 
excitations. It is worth repeating here that physically A encodes information 
about properties of the series (fT^ and has the same theoretical status as the 
corresponding np condensates in the following sense: knowledge of all terms 
in (1141) would allow to reconstruct A. Of course, it is impossible in practice 
and one has to use other methods to study large distance asymptotics of 
FC. On more phenomenologic level it was discussed in ^6], while rigorous 
definition was given in the framework of the FC methoco (see [7], [27] and 
[TT] for review). 

It is important to find also the world-line representation for gluon Green's 
function, which is done in the framework of background perturbation theory 
(see [28] for details). In this way one finds (where a, b are adjoint color indices 
and dash sign denotes adjoint representation) 

ab 



/ ds{Dz)^ye-''PaeMw A^dz^)P^{x,y, s)\ , (16) 

-'0 Jy ) 

Pj:{x,y,s) = PFexp{2ig / Fxa{z{T))dT). 

Jo 



where 



^ Often referred to as Stochastic Vacuum Model 
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All the above reasoning about asymptotic expansions of the corresponding 
W-loops is valid here as well. 

Thus the central role in the discussed method is played by quadratic 
(Gaussian) FC of the form 

D^^^M = 9\TTF,,{xMx,y)F,^{yMy,x)), (17) 

where F^^iy is the field strength and $(x, y) is the parallel transporter. Corre- 
lation lengths Aj for different channels are defined in terms of asymptotics of 
f lT7|) at large distances: exp (— |x — y\/Xi). The physical role of Aj is very im- 
portant since it distinguishes two regimes: one expects validity of potential- 
type approach describing the structure of hadrons of spatial size R and at 
temporal scale Tg for Aj <^ R, Tg, while in the opposite case, when Aj ^ R, Tg 
the description in terms of spatially homogeneous condensates can be applied. 

At zero temperature the 0(4) invariance of Euclidean space-time hold^l 
and FC fll7p is represented through two scalar functions D{z), Di{z) (where 
z = X — y) as follows 



d 

va 



dz^ 



D,{z). (18) 



One has to distinguish from the very beginning perturbative and np parts 
of the correlators D{z), Di{z). Beginning with the former, one easily finds 
at tree level 

D^'\z) = ; Dfiz)=C,if)^\ (19) 



where C2(/) is fundamental Casimir C2{f) = {N^ — l)/2Nc. At higher orders 
situation becomes more complicated. Namely, one has at n-loop order 

DP'"(z) = Df{z) ■ G(")(^) ; ^r(^) = Df{z) ■ G'('\z) (20) 

where the gauge-invariant functions G'^^'^z), G^^^z) has the following general 
structure: 

G(")(^),G(")(z)~a: [c„(log/iz)" + ...] 



^AU treatment in this paper as well as averaging over vacuum fields is done in the 
Euclidean space. Notice that only after all Green's functions are computed, analytic 
continuation to Minkowskii space-time can be accomplished. 
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where c„ is numerical coefficient, is renormalization scale and we have 
omitted subleading logarithms and constant terms in the r.h.s. Explicit ex- 
pressions for the case n = 1 can be found in |29l |30] . 

Naively one could take perturbative functions DP'"'{z), Di'^{z) written 
above and use them for computations of W-loop or static potentials in Gaus- 
sian approximation. This however would be incorrect. The reason is that 
at any given order in perturbation expansion over one should take into 
account perturbative terms of the given order coming from all PCs, and not 
only from the Gaussian one. For example, at one loop level the function 
DP'^{z) is nonzero which naively would correspond to area law at perturba- 
tive level. Certainly this cannot be the case, since self-consistent renormal- 
ization program for W-loops is known not to admit any terms of this sort 
[31]. Technically in FC language the correct result is restored by cancela- 
tion of contributions proportional to D^'^(z) to all observables by the terms 
coming from triple FC (see details in |32j). 

In view of this general property of cluster expansion it is more natural 
just to take relation D^'^{z) = as valid at arbitrary n, having in mind that 
proper number of terms from f[T^ has to be included to get this result. Thus 
we assume that the following decomposition takes place: 

D{z) = D-^iz) ; D,iz) = D{{z) + D7{z) (21) 

and D{z) has a smooth limit when 2; — * 0. As for Di{z) its general asymptotic 
at 2 — 5> reads: 

D,{z) = -^ + % + 0{z') (22) 

Z^ Z'^ 

where c and 02 weakly (logarithmically) depend on z. Notice that by -Di(O) 
we always understand np "condensate" part in what follows: 

G2 = ^(D"^(0) + 1^7(0)) (23) 

The term l/z^ deserves special consideration. One may argue, that at 
large distances the simple additivity of perturbative and np contributions to 
Vqq potential is violated [33]. It is interesting whether this additivity holds 
at small distances and, in particular, mixed condensate 02 was suggested in 
[31] and argued to be phenomenologically desirable. 

As we demonstrate below our approach is self-consistent with 02 = 0, 
i.e. when perturbative- additivity holds at small distances and there is 
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no dimension-two condensate. However, at intermediate distances a compli- 
cated interrelation occurs, which does not exclude fl22l) being relevant in this 
regime. 

We proceed with general analysis of f|T7j) . As is proved in [7J, D{z), Di{z) 
do not depend on relative orientation of the plane (/iz/), plane (Act), and the 
vector Za- This orientation can be of the following 3 types: a) planes (fiu) 
and (Act) are perpendicular to the vector z^, b) planes (/xz^), (Act) are parallel 
or intersecting along one direction and Za lies in one of planes; c) planes {^u) 
and (Act) are perpendicular and Za lies necessarily in one of them. 

It is easy to understand that this classification is 0(4) invariant, and 
therefore one can assign indices a, b, c to D^,^ xa, resulting in general case in 
three different functions. Physically speaking the case a) refers to e.g. color 
magnetic fields Fik, Fim with Za in the 4th direction, and the corresponding 
FCs will be denoted as D±{z). The case b) refers to e.g. the color electric 
fields Ei{x), Ek{y), connected by the same temporal links along 4th axis, and 
the corresponding FCs are denoted as D^\{z). Finally, in case c) only Di part 
survives in f[T5]) and it will be given by the subscript EH, Df^. 

In case of zero temperature, when 0(4) invariance holds, these three 
functions D±{z), D\\{z), Df^ (z) are easily expressed via D{z), Di{z) 

D^{z) = D{z) + Di{z), 
DJz) = D{z) + Di{z) + z- 



,dD,iz) 



dz^ 



Df"{z) ^ D,{z), 



(note that just D±,D^\ were measured on the lattice in [14]). For nonzero 
temperature the correlators of colorelectric and colormagnetic fields can be 
different and in general there are five FCs: D^{z), Df{z), D^{z), D^{z), and 
Df^[z). As a result one obtains the full set of five independent quadratic 
FCs which can be defined as follows 

/(M,(x)$if,(y)$+) = 6,,{D^{z) + D^{z) + z^^) - (24) 



g'{TTE,ix)<!>E,iy)^+) = 5,,{D^{z) + (z) + zl^) + z,z,^, (25) 



QJjEH 

g\i:iHi{x)^E,{y)<^+) = djkZkZ^^^. (26) 
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It is interesting that the structure (12^ - fl26p survives for nonzero tempera- 
ture, where 0(4) invariance is violated: it tells that the equations 
give the most general form of the FCs. Note that at zero temperature, 
= D^,D^ = Df at coinciding point; the 0(4) invariance requires that 
colorelectric and colormagnetic condensates coincide: g"^ {Tr Hi{x)Hi{x)) = 
g'^{TTEj{x)Ej{x)), hence 

Z}^(0) + Df(0) = D^(0) + Df(0). (27) 

Here we assume the np part of all FCs is finite. 

At 2 7^ for zero temperature FCs can be expressed through Dj_{z), [z) 
which do not depend on whether Ei or Hk enters in them, since e.g. can 
be transformed into by an action of 0(4) group elements (the same for 
Df, D^). From this one can deduce that both Di and D do not depend on 
subscripts E, H for zero temperature and Df^ coincides with Di. 

Below we shall illustrate this coincidence by concrete calculations of 
jjE ^ jjH through the gluelump Green's functions and show that the cor- 
responding correlation lengths satisfy the relations (with obvious notations): 

A^ = A^ = A; Af = Af = Af^ = Ai. (28) 

As is shown below, all correlation lengths A^ appear to be just inverse masses 
of the corresponding gluelumps Xj = 1/ Mj . 

Having analytic expressions for FCs one might ask how to check them 
versus experimental and lattice data. On experimental side in hadron spec- 
troscopy one measures masses and transition matrix elements, which are 
defined by dynamical equation and the latter can be used of potential type 
due to smallness of Xj. In static potential only integrals over distance enter 
and the spin-independent static potential can be written as [3 153] 

pr POO 

VQQ{r) =2 (r - X)dX / duD^iVX^ + + 
Jo Jo 

pr POO 

+ XdX Df(VA2 + z/2). (29) 
Jo Jo 

At this point one should define how perturbative and np contributions 
combine in D^,Df, and this analysis was done in [32]. Making use of f[T9l) 
together with definition of the string tension [7], 

= ^ 1 J d'zD^iz), (30) 
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one obtains from fl29l) the standard form of the static potential for = 3 at 
distances r 

VQ^{r) = a\-^-^ + 0{X^/r,al). (31) 

As is argued below, ~ 0.1 Fm, so that the form flHT]) is applicable in the 
whole range of distances r > 0.1 fm provided asymptotic freedom is taken 
into account in as{r) in (13T|) . 

At the smallest values of r, r < A^, one would obtain a softening of 
confining term, err — > cr^ [7], which is not seen in accurate lattice data at 
r > 0.2 fm [Ml ETj, imposing a stringent limit on the value of A^; A^ < 0.1 
Fm, see discussion in [37] . 



3 Spin-Dependent Potentials and FC 

We now turn to the spin-dependent interactions to demonstrate that FC can 
be extracted from them [351 HHl HI] . Wl'^)) ^2(^)) ^3(^)5 ^4(^^)5 plus static 
term V"Qg(r) contain in integral form all five FC: -D^, Df, , , D^^ and 
one can extract the properties of the latter from the spin-dependent poten- 
tials. This procedure is used recently for comparison of analytic predictions 
[20] with the lattice data [19]. 

To express spin-dependent potentials in terms of FC, one can start with 
the W-loop ([8|), entering as a kernel in the qq Green's function, where quarks 
q, q can be light or heavy. Since both kernels Pa, contain the matrix 
o'fiuF'^u, the terms of spin-spin and spin-orbit types appear. As before, we 
keep Gaussian approximation, and in this case spin- dependent potentials can 
be computed not only for heavy, but also for light quarks. They have the 
Eichten-Feinberg form [32]: 



ridiag)fjy^ _ ( ^iLi _ o^lA f ]_de_ 2dVi{R) 
Afij A^l j \RdR RdR 



^ a2Li - aiL2 I dV2 ^ aia2V4(^) ^ (3ai^a2^ - o'io'2^^)^3 ^^2) 
2yUi/i2 R dR 12/ii/i2 12/ii/i2-R2 
Spin-spin interaction appear in W^, Eq.® which can be written as 

L U«E a«B U(^)E a(2)B 
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where Zi^2 = z{ti^2) and spin-orbit interaction arises in fll7l) from the products 

It is clear, that the resuhing interaction will be of matrix form (2 x 2) x 
(2 X 2) (not accounting for Pauli matrices). If one keeps only diagonal terms 
in a^yF^jj (as the leading terms for large /i, ~ M) one can write for the 
spin-dependent potentials the representation of the Eichten-Feinberg form 



At this point one should note that the term with ^ in (!32|) was obtained 
from the diagonal part of the matrix {m — D)a^yFpu, namely as product 



idkDk ■ (TiEi, see [35J HQ] for details of derivation, while all other potentials 
Vi,i = 1,2,3,4 are proportional to FCs {Hi^Hk^) . One can relate PCs of 
colorelectric and colormagnetic fields to D^, Df, , defined by fl24 
fl26|) with the following result: 



1 dVi 
R~dR 



■oo Jo ^ \ 



1^ 

R dR 



dv 



MX 



H 



duR 



2dDi'{R,u) 



dR^ 



V, = J dz./3D^(i?,z/) + 3Df(i?,z/) + 2i?2_^J 



(33) 

(34) 
(35) 
(36) 



1 d6{R) 
R dR 



^ 



D^{X,u) + Df{X,u) + {X' + u')^ 



One can check, that the Gromes relation [18] acquires the form 

V;{R)+e'{R)-V^{R) = 

/oo r pR -| 

du / dA(D^(A,z/)-D^(A,z/)) + -i?(Z}f(i?)-Z}f(i?)) 
-oo L"'0 2 



(37) 



(38) 



Por T = 0, when = , Df = , the Gromes relations are satisfied 
identically, however for T > electric and magnetic correlators are certainly 
different and Gromes relation is violated, as one could tell beforehand, since 
for T > the Euclidean 0(4) invariance is violated. 
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4 Gluelumps and FC 



In this section we establish a connection of FC D{z) and Di{z) with the 
gluelump Green's functions. 

To proceed one can use the background field formalism [28] [13], where 
the notions of valence gluon field and background field 5^ are introduced, 
so that total gluonic field is written as 

A^ = af, + B^. (39) 

The main idea we are going to adopt here is suggested in the second paper 
of [15]. To be self-contained, let us explain it here in simple form. First, we 
single out some color index a and fix it at a given number. Then in color 
components we have, by definition, 

= + (1 - ns; 

and there is neither summation over a in the first term nor over c in the 
second. As a result the integration measure factorizes: 

c Cy^a 

SO one first averages over the fields B'^, and after that one integrates over 
the fields a^. The essential point is that the integration over 'DB'^ provides 
colorless adjoint string attached to the gluon a^, which keeps the color index 
"a" unchanged in the course of propagation of gluon a in background made 
of B. This is another way of saying that gluon field of each particular color 
enters QCD Lagrangian quadratically. Despite such separation of the mea- 
sure is a kind of trick, the formation of adjoint string is the basic physical 
mechanism behind this background technic and it is related to the properties 
of gluon ensemble: even for = 3 one has one color degree of freedom 
and 7 fields S^, also confining string is a colorless object and therefore the 
singled-out adjoint index "a" can be preserved during interaction process of 
the valence gluon with the background (-B^). These remarks to be used in 
what follows make explicit the notions of the valence gluon and background 
field. 

We assume the background Feynman gauge D^a^ = and assign field 
transformations as follows 

a^^U+a^U, B^^U+{B^ + -d^)U. (40) 
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As a result the parallel transporter y) = Pexpig B^{z)dz keeps its 
transformation property and every insertion of between $ transforms 
gauge covariantly. In what follows we shall assume that $ is made entirely 
of the field B^. Note that in the original cluster expansion of the W-loop 
W{B + a) in FC [13]. 

W{B + a) = (TrP exp tg I {B^ + a^)rfz^) (41) 

Jc 

it does not matter whether {B^ + a^) or B^ enters ^'s since those factors 
cancel in the sum. Thus we assume from the beginning that $'s are not 
renormalized, since such renormalization is unphysical anyway according to 
what have been said above. 

As for renormalization having physical meaning, it is known [31] to reduce 
to the charge renormalization and renormalization of perimeter divergences 
on the contour C. For static quarks the latter reduces to the mass renormal- 
ization. 

In background field formalism one has an important simplification: the 
combination gB^ is renorminvariant [28l [31] , so is any expression made of 
Bfj^ only. This point is important for comparison of FC with the lattice 
correlators in [TSl ED] • Renormalization constants for field strengths Zf, used 
there account for finite size of plaquettes and they are similar to the lattice 
tadpole terms. 

Now we turn to the analytic calculation of FC in terms of gluelump 
Green's function. To this end we insert f l39l) into f[T71) and have for F^i,{x): 

F^y{x) = d^A^ - dyA^ - ig[Af,, Ay] = D^^a^ - D^a^ - ig[af„ a^] + F^f \ (42) 

Here, the term F^f'' contains only the field B'^. It is clear that when one 
averages over field a° and sums finally over all color indices a, one actually 

( Fi) ( R) 

exploits all the fields with color indices from F^^y , so that the term F/iy can 
be omitted, if summing over all indices a is presumed to be done at the end 
of calculation. 

As a result D^y \„ can be written as 

D,y,x. {x,y)= D + « A. + (43) 

where the superscript 0,1,2 refers to the power of g coming from the term 
ig[a^,ay]. 
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We can address now an important question about the relation between 
FCs and gluelump Green's functions. We begin with 1-gluon gluelump, whose 
Green's function reads 

= (Tr,a^(x)l>(x,i/)a,(y)). (44) 

According to fj40l) G^u? ^ (x, y) is a gauge invariant function. 

As shown in |T5j, the first term in (jl3D is connected to the functions 
, and it can be written as follows 

D^M^^.y) = ^ {a^^^- (^'^) +^™} + <A., (45) 

where ^^^Ixa contains contribution of higher FCs, which we systematically 
discard. 

On the other hand one can find G^^u^ {x, y) from the expression |21j written 

as 

poo 

G%^'\x,y) = TTj ds{DzUeM-K){W;:AC.y)). (46) 
where the spin-dependent W-loop is 

K^{C,y) = PPf |exp(z(7 j Bxdzx) exp f| (47) 



and the gluon spin factor is expF = exp(2z(7 drFBiz^r))) with Fb made 
of the background field only. 

Analogous expression can be constructed for Green's function of 2-gluon 
gluelump. It is given by the following expression 

G%%ix,y) = {TT,{r'j''^^a';{x)al{x)nix,y)T^x 

X aiiy)aliy)) ^ Nl{Nl - l){5,y5,, - 5,J,^)G^^^'\z). (48) 
At small distances G'*^^^'-*(x, y) is dominated by perturbative expansion terms, 

G^'''\z) ~ 1, (49) 

however, as we already discussed in details, all these perturbative terms are 
canceled by those from higher FCs (triple, quartic, etc.), therefore expansion 
in fact starts with np terms of dimension four. 
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To identify the np contribution to G^'^^^\z) we rewrite it as follows: 

roQ poo 

G^^3'\z)= dsj ds2iDz,)oADz2)o.TrW^{C^,C2) (50) 
Jo Jo 

where the two-gluon gluelump W-loop Wy:{Ci,C2) is depicted in Fig. 6 and 
can be written as (in the Gaussian approximation) 

Ws(Ci,C2) = exp|-i^y" d7i^,iu)d7r,^{v){F^,iu)^Fx,iv))y (51) 

and the total surface S consists of 3 pieces, as shown in Fig. 7 

F^^diT^u{u) = F^yds^u{u) - 2igdT{F{u)) (52) 

where (F) has Lorentz tensor and adjoint color indices F°'l' and lives on gluon 
trajectories, i.e. on the boundaries of Si. In full analogy with (145!) we have 

D{z) = ^'(^^-^) g(2.0(^) (53) 

The crucial point is that the Green's function (H5|) . (150|) can be calculated 
in terms of the same FCs D{z), Di{z). Indeed one has 

G'^-\z) = {f\expi-H^\z\)\z) (54) 

where the index w stays for 1-gluon or 2-gluon gluelump Hamiltonians. The 
latter are expressed via the same FCs D{z), Di{z) (see [45] and references 
therein) : 

= Ho [fj] + AHl [/i, u] + AHcoui [Di] + AHst„ng [D, u] (55) 
where the last term Hstring[F>, v\ depends on Di^z) via adjoint string tension 

= - \ d^zD{z) (56) 

and Hamiltonian depends on einbein fields fi and u. The term corresponding 
to perimeter Coulomb-like interaction AHcoui[Di] depends on Di{z) (com- 
pare with (1291). (|32|)). 

The self-consistent regimes correspond to different asymptotics of the 
solutions to these equations. In Coulomb phase of a gauge theory both Hig, 
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exhibit no mass gap, i.e. large z asymptotic of fl54l) is power-like. The 
function D{z) vanishes in this phase and W-loop obeys perimeter law. In 
the confinement phase realized in Yang-Mills theory at low temperatures a 
typical large-z pattern is given by 

D{z),Dl{z)^eM-\AI^^) 

i.e. there is a mass gap for both Hig, H2g. 

Confining solutions are characterized by W-loops obeying area law. In 
other words, Hamiltonians expressed in terms of interaction kernels depend- 
ing on D(z), Di{z) exhibit mass gap if these kernels are confining. On the 
other hand the same mass gap plays a role of inverse correlation length of 
the vacuum. This should be compared with well known mean-field technique. 
In our case the role of mean- field is played by quadratic FC, which develops 
non-trivial Gaussian term D{z). Notice that the exponential form of its large 
distance asymptotics exp(— |2;|/A) (and not, let's say, exp^—z^/X^) ) is dic- 
tated by spectral expansion of the corresponding Green's function at large 
distances. 

Full solution of the above equations is a formidable task not addressed 
by us here. Instead, as a necessary prerequisite we check below different 
asymptotic regimes and demonstrate self-consistency of the whole picture. 

We begin with small distance region. The np part of contribution to 
Di{z) at small distances comes from two possible sources: the area law term 
(first exponent and the expF term in (H7|) ). Both contributions are depicted 
in Fig.l and Fig.2 respectively. We shall disregard the term expF in this 
case, since for the one-gluon gluelump it does not produce hyperfine inter- 
action and only gives rise to the np shift of the gluon mass, which anyhow 
is eliminated by the renormalization (see Appendix for more detail). This is 
in contrast to the two-gluon gluelump Green's function generating D^, , 
where the hyperfine interaction between two gluons is dominating at small 
distances. 

As a result Eq. without the (exp F) term yields 



A(.) = ^^^ + -G,. (57) 



It is remarkable that the sign of the np correction is positive. 

At large distances one can use the gluelump Hamiltonian for one-gluon 
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gluelump from [17] to derive the asymptotics 

D,(.) = ^^^^^^^;j^°'^^"^^- e-<^N, |.|Mr»l. (58) 

where M^'^^ = (1.2 ^ 1.4) GeV for af = 0.18 GeV^ 

We now turn to the FC D{z) as was studied in JW\. The relation (1531) 
connects D{z) to the two-gluon gluelump Green's function, studied in [47] 
at large distances. Here we need its small - z behaviour and we shall write 
it in the form 

G^^sO (^) = ^0 (z) + Ggf' ) (z) (59) 

where G^^^\z) contains purely perturbative contributions which are sub- 
tracted by higher-order FCs, while Gnp^\z) contains np and possible pertur- 
bative - np interference terms. We are interested in the contribution of the 
FC {FF) to G^'^^^\z), when z tends to zero. 

One can envisage three types of contributions: 

a) due to product of surface elements ds^yds\„, which gives for small 

~^I/v" ) ' similar to the situation discussed for Di. 
This is depicted in Fig. 7. 

b) contribution of the type dTidT2{FF) , where dr and dT2 belong to dif- 
ferent gluon trajectories, 1 and 2, as depicted in Fig.8. This is the hyper- 
fine gluon-gluon interaction, taken into account in [l^ in the course of the 
gluelump mass calculations (regime of large z), however in that case mostly 
the perturbative part of (FF) contributes (due to Di). Here we keep in (FF) 
only the NP part and consider the case of small z. 

c) Again the term dr^drj, but now i = j. This is actually a part of the 
gluon selfenergy correction, which should be renormalized to zero, when all 
(divergent) perturbative contributions are added. As we argued in Appendix 
of [16], we disregard this contribution here, as well as in the case of Di{z) 
(one-gluon gluelump case). Now we treat both contributions a) and b). 

a) In line with the treatment of Di, Eq. (!57|) . one can write the term, 
representing two-gluon gluelump as two nearby one-gluon gluelumps which 
yields for D{z), 

AD^iz) = (60) 

b) in this case one should consider the diagram given in Fig. 8, which yields 
the answer (for details see Appendix). 
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^(2.0 . . ^ r d^wd^w'D{w-w') 

^ ^ ' N^-lJ {A7T^y{w-z)V{w'-z)V 

which contributes to D{z) as 

ADb{z) = 2Nlg%z). (62) 

where at small z, h{z) ~ log^ (^T^)' order of correlation 

length A. 

At large distances one use the two-gluon gluelump Hamiltonian as in 
[18] and find the corresponding spectrum and wave functions, see [IE] and 
appendix 5 of [16] for details. As a result one obtains in this approximation 

D{z) = ^'^^l~^\ .la)e-<"'\^\ Mf|z|»l (63) 

where Mq is the lowest two-gluon gluelump mass found in [17] to be about 
Mq^^^ = (2.5 ^ 2.6) GeV 

We shall discuss the resulting properties of D{x) and -Di(x) in the next 
section. 



5 Discussion of Consistency 

We start with the check consistency for D{z). As is shown above, D{z) has 
the following behavior at small z 

D{z) ^ -m,aMz))G, + N^^^^D{Xo) log' (^) (64) 

Since as{fi{z)) ~ 27r//?o log(A2;)~^ the first term is subleading at 2; — and 
the last term on the r.h.s. tends to a constant: 

Since from fl65l) one can infer, that -D(O) ~ 0.15D{Xq) for Nc = 3, where 
-^0 ~ -^^ ■ So D{z) is an increasing function of z at small z, z < Aq, and for 
z ^ X one observes exponential falloff. The qualitative picture illustrating 
this solution for D{z) is shown in Fig. 5. 
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This pattern may solve qualitatively the contradiction between the values 
of -D(O) estimated from the string tension D„{0) ^ ^ ^ 0.35 GeV^ and the 
value obtained in naive way from the gluon condensate DgjIO) = fs^s ~ 
(0.007H-0.012) GeV^. One can see that D^{0) ^ {30^54:)DgM. This seems 
to be a reasonable explanation of the mismatch discussed in the introduction. 

As was shown in |T6] , the large distance exponential behavior is selfconsis- 
tent, since (assuming that it persists for all z, while small z region contributes 
very little) from the equality a = TTX'^Dcr{0), comparing with (|63l) one obtains 

0.1 . 87:W,{N^ - = (66) 

where in Q;s(/i) the scale /i corresponds roughly to the gluelump average 
momentum (inverse radius) /iq ~ 1 GeV. Thus ( 166!) yields as{f^o) ~ 0.4 
which is in reasonable agreement with «<, from other systems 
We end this section with discussion of three points: 

1) D{z) and Di{z) have been obtained here in the leading approximation, 
when gluelumps of minimal number of gluons contribute: 2 for D{z) and 1 
for Di{z). In the higher orders of 0{as) one has an expansion of the type 

D{z) = D^^^'\z) + cialD^^'^'\z) + C2alD^^''\z) + ... 

D^{z) = D^^'J'\z) + d^alD^^''\z) + ... (67) 

Hence the asymptotic behavior for D{z) will contain exponent of Mq^''|2| too, 
but with a small preexponent coefficient. 

2) The behavior of D[z), D^^^\z) at small z is defined by NP terms of 
dimension four, which are condensate G2 as in f|T2l) . and the similar term 
from the expansion of exp F, namely, 

(exp F) = 1 + 4// dr dr' {F{u{t))F{u{t')) + ... (68) 
Jo Jo 

therefore one does not encounter mixed terms like O(^), however if one 
assumes that expansion of {F{0)F{z)) starts with terms of this sort as it 
was suggested in |3l|, then one will have a selfconsistent condition for the 
coefficient in front of this term. 

3) To study the difference between , at T < Tc one should look 
at Eq. fH5l) and compare the situation, when /i = A = 4, u = a = i 
and take z = x — y along the 4-th axis (for -Df ), while for one takes 
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H = X = i, i' = a = k and the same z. One can see that in both cases one 
ends up with the one-gluon gluelump Green's function G^,^ fll^ which in the 
lowest (Gaussian) approximation is the same Gf^^ = Sfiuf{z), and f{z) is the 
standard lowest mass (and lowest angular momentum) Green's function. 

Hence Df = in this approximation and for T = this is an exact 
relation, as was discussed above. Therefore 

Ai = Af = Af = 1/M(5^\ M(5^^ ^ 1.2 ^ 1.4 GeV, Ai ^ 0.2 ^ 0.15 fm. (69) 

The value M^'^^ in 0691) is taken from the calculations in |47j . The same is true 
for D^,D^, as it is seen from fHS]) . where G*^^^') corresponds to two-gluon 
subsystem angular momentum L = independently of /iz/, Act. 
Hence one obtains 

A = A-^ = A^ = A = 0.08 fm, M^^^ ^ 2.5 GeV. (70) 

where the value Mq ^ 2.5 GeV is taken from [47] . 

6 Conclusions 

We have derived, following the method of [28] the expressions for FC D{z), 
Di{z) in terms of gluelump Green's functions. This is done in Gaussian 
approximation. The latter are calculated using Hamiltonian where np dy- 
namics is given by D"'P{z), D^^{z). In this way one obtains selfcoupled equa- 
tions for these functions, which allow two types of solutions: 1) D'"'P{z) = 0, 
Di^{z) = 0, i.e. no np effects at all; 2) D"'P{z), D^^{z) are nonzero and 
defined by the only scale, which should be given in QCD, e.g. string tension 
a or Aqcd- All other quantities are defined in terms of these basic ones. 
We have checked consistency of selfcoupled equations at large and small dis- 
tances and found that to the order 0{as) no mixed perturbative-np terms 
appear. The function Di{z) can be represented as a sum of perturbative and 
np terms, while D{z) contains only np contributions. 

We have found a possible way to explain the discrepancy between the 
average values of field strength taken from G2 and from a by showing that 
D{z) has a local minimum at z = and grows at 2; ~ A. Small value of A and 
large value of the gluelump mass Mgi = 1/A ~ 2.5 GeV explains the lattice 
data for A ~ 0.1 Fm. Thus the present paper argues that relevant degrees of 
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freedom ensuring confinement are gluelumps, described self-consistently in 
the language of FCs. 
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A Nonperturbative contributions to Di(z) and 
D{z) at small z 

We start with Di{z) expressed via gluelump Green's function flUl) . (jlH]) 

The leading NP behavior at small z, proportional to {F{x)F{y)), is ob- 
tained from the amplitudes, shown in Fig.l and Fig.2. For the amplitude 
of Fig.l one can use the vacuum average of the W-loop (fT2|) in the limit of 
small contours and neglecting the factor expF in ( P7|) . One obtains for the 
adjoint loop 

(WiC^y)) = exp (-^^^^^7^^) (71) 

where 7 = ^C2U) ~ n^~i (o^igi^^^'l derivation in [16] referred to the fun- 
damental loop). Proceeding as in Appendix 2 of p!6] one obtains the NP 
correction ADi 

2 2 

AD,a = C,{ad3)^G, = ^-G^ (72) 

where we have introduced the standard definition G2 = ^{F%iS^)F^uiS^)) ■ 
One can check, that ADi = C(A°). 

We turn now to the amplitude of the Fig.2. The corresponding f{x) for 
a; — >■ can be written as follows 

A/ib(a;) = J d\d\GiO,u)G{u,v)G{v,x)A7i^N^G2 (73) 



where the gluon Green's function G{x,y) for small \x — y\ can be replaced 
by the perturbative part: G{x,y) Go{x,y) = i^^l_yy2 - 
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Taking into account that ADi{x) = — ^^A/(x), one obtains for x —>■ 

0, 

ADM = ^g'G^Hx) (74) 

where we have defined 

T(x) = _ A f ^'^^'^ (75) 

^ ' dx^ J {ATi^Yu'^{u-v)\v-xY' ^ ' 

The integral /(O) diverges at small v and large m, v. The latter divergence 
is removed since at large arguments G{x,y) is damped by confining force, 
keeping the propagating gluon nearby the 4-th axis, x^ = = where the 
static gluonic source resides. One can estimate I{x ^ 0) ~ 0((/n^)")). 
Since G2 is connected to D, Di 

G2 = -^(D^iO) + DfiO) + D^(0) + Df (0)) (76) 

one can see in fl74p . that coefficients of Df on both sides of fl74p have the same 
order of magnitude and the same sign, suggesting a selfconsistency on this 
preliminary level. However, as we argued before and in [16] this contribution 
is actually gluon mass renormalization, which is zero. This is especially clear 
when (FF) in Fig. 2 is replaced by its perturbative part. 

We turn now to the most important case of the confining FC D^{x). The 
corresponding two-gluon gluelump amplitude is depicted in Fig.3, and the 
np contributions to D^{x) are given in Fig.4 and Fig.5. 

We start with the amplitude of Fig.4, which can be represented as a 
doubled diagram of Fig.l and therefore the contribution to the two-gluon 
gluelump Green's function G^'^^^\x) will be 

G^^3i)u) = \ 1^ (77) 

^' {A-K^x^f 24iVe(47r2)- ^ ' 

Now using relation f[T21) one obtains for the np contribution of Fig.4 

AD.(x) = ^-^^^AG^^^Hx) = (78) 

Note, that this contribution is finite at a; — > 0, however with the negative 
sign. We now turn to the np contribution of Fig.5, which can be obtained 
from fH5]) inserting there the product of operators Ffg{w)d^wFf/gi{w')d'^w' , 
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where fg{f'g') are adjoint color indices, Ffg = F°-Tjg = F°-{—i)fafg- Using 
the formulas 

fabcfa'bc = N^5aar, f abc fade hdf fee, = -fSfg, (79) 

one arrives at the expression 

AD,{x)=2N^g^h{x), (80) 

h{x) = I d^wd'^w'Diw - w')G{w)G{w')G{w - x)G{w' - x) 



where G{y) is the gluon Green's function in the two-gluon gluelump; the 
gluon is connected at large distances by two strings to another gluon and to 
the static gluon source, see Fig.5. At small x (we take it along axis 4 for 
convenience) the integral in h{x) grows when G{y) is close to the 4th axis 
and becomes the free gluon, G{y) Go{y) = ^^2^2 • As a result one obtains 

f d^wd'^w' D(w — w') , . 

h[x) ho{x) = / r^— 81 

and one should have in mind that the integral for x — diverges both at 
small and at large w, w' . The small w, w' region will give the terms 0{ln-%), 
and at large w, w' the integral is protected by the fall-off of G{y) at large 
y due to confinement, therefore we must imply in the integrals in ( 1811) the 
upper limits for w,w' at some Aq ^ , ^ . Then introducing for w(w, W4) 

polar coordinates |w| = p sin 6,W4 = p cos 6 (and the same for w'), one arrives 
at the integral 

r sin^^d^ _ ^ { \ 
Jo p^-2xpcos9 + x^~2^^P-''^ ^^^^ 

or (2; > p), and finally one has the estimate of contribution to fISUl) from 
the region of small w, w'{w\w\, \w'\ < Aq) 

^^^"^"4(4^^°^^ J ^ ^ 

and finally for D{x) one obtains (omitting the perturbative term in (1771) . 
which cancels. 



BW = -#G. + ^a2fl(A.)V(^|. (84) 
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Note, that the last term on the r.h.s. dominates at small x and compen- 
sates the decrease of the a1 term, since 



X 



Hence 

D(a:<Ao)^^(^) D{\o)<D{\o) (86) 
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